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Influence of a Ridge on the Low-Frequency 

Ground Wave 

James R. Wait and Anabeth Murphy 

The problem of a plane wave incident on a semielliptical boss on an otherwise perfectly 
conducting flat ground plane is considered. A solution in terms of elliptic wave functions is 
obtained. Numerical values of the field on the near and far side of this idealized ridge are 
given for a base width of about two-thirds of a wavelength and various ellipticity ratios. 

It is now relatively straight forward to make predictions of ground wave propagation on a 
smooth, homogeneous, spherical earth. 1,2 Some progress has also been made recently in de- 
vising techniques to calculate the fields over an inhomogeneous and nonsmooth earth. For 
example, at high frequencies, the effect of mountains or ridges obstructing the transmission 
path can be treated by methods of physical optics. 3 Another approach, disclosed very re- 
cently, simulates the obstacle by four spherical surfaces. 4 The transmission loss over each 
segment is then added, being an approximation valid at very high frequencies in most cases. 

At low frequencies (less than 1 Mc) obstacles such as ridges are, no longer, large compared 
to the wavelength, and it is necessary to use a different approach. A suggested model to study 
the effect of a nonsmooth ground for low frequencies is a semielliptical boss on an otherwise 
flat ground plane. Assuming that the source of vertically polarized waves is at a distance from 
the elliptic cylinder, large compared to the wavelength, it is sufficient to consider the incident 
wave to be plane. Furthermore, at low frequencies the attenuation of the incident wave is 
negligible for a distance of several wavelengths on either side of the ridge. For this reason, the 
ground plane and the surface of the ridge is considered to be perfectly conducting. 

The semielliptical perfectly conducting boss is shown in figure 1. With reference to a carte- 
sian coordinate system (;x, y, z) the ground plane is y=0 and the surface of the boss is x 2 /b 2 + 
y 2 /a 2 =l y where b and a are the semimajor and semiminor axes of the ellipse, respectively. 
Elliptic-cylinder coordinates (u, v, z), which are confocal with the elliptical boss, are connected 
to the cartesian coordinates by 

x=c cosh u cos v 



y=c sinh u sin v 

Z=Z, 



(i) 



where c, the semifocal distance, is given by 

c=(b 2 -a 2 ) y \ 
The magnetic field of the incident wave has only a z component and is given by 

where H is the amplitude and /3=27r/wavelength. Employing an addition theorem, 5 this can 
be rewritten in terms of a particular solution of the wave equation in elliptic-cylinder coordi- 
nates as follows* 

//'- =# (8t)' i S %£$ Je m {c,u)e imr '\ (2) 



' K. A. Norton, "The calculation of the ground wave field intensity over a finitely conducting spherical earth, Proc. Inst. Radio Engrs. ?9, 
623 (December 1951). 

2 U. Bremmer, Terrestial radio waves (Elsevier Publishing Co., Ltd., New York, N. Y., 1949). 

3 F. H. Dickson, J. J. Egli, J. W. Herbstreit, and G. S. Wickizer, Large reductions of V. H. F. transmission loss and fading by the presence 
of a mountain obstacle in beyond line of sight paths, Proc. Inst. Radio Engrs. 41, 967 (August 1953). See also a discussion of this paper by J. H. 
Crysdale, Proc. Inst. Radio Engrs. 43, 627 (May 1955). 

4 K. A. Norton, P. L. Rice, and L. E. Vogler, The use of angular distance in estimating transmission loss, Proc. Inst. Radio Engrs. 43, 1488 
(October 1955). 

s P. M. Morse and H. Feshbach, Methods of theoretical physics 2, 1572 (McGraw Hill Book Co., New York, N. Y., 1953). 
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Figure 1. Semielliptical boss. 

where Se m is the even angular Mathieu function, Je m is the even radial Mathieu function of 
the first type, and Nm is the normalization constant for Mathieu functions of order m. 

The secondary or scattered field, H sec ', is now chosen to satisfy the boundary condition that 
the normal derivative of the total field is zero at u=Ui, the surface of the elliptic cylinder, and 
at v = and w, the surface of the ground plane. Furthermore, the secondary field must contain 
the radial Mathieu function He® } of the fourth type to comply with the radiation condition at 
infinity. It then follows that 



H™- = -H (8w)KTi 



Sejc,v) 



=o N e m (c) 



' i)Je m (c,u) 
du 

dHe™(c,u) 
du 



He^(c,u)e imT/2 . 



(3) 



The Mathieu functions are not tabulated directly but must be calculated from series expansions 
whose coefficients are available in tabular form. 6 In the present notation 



and 



Se 2 r + p(c,v)=J2, De 2k + P cos [(2k+p)v] 

fc=0 



Je»+ P (fi,u)=( -) (— l) r ]C (—l) k De 2 j c + P J2k+p(c cosh u), 



(4a) 



(4b) 



where r is an integer, p = Q or 1, J 21c + P is the Bessel function of the first type or order 2k+p, and 
De 2k+P are coefficients that are functions of 2r+p and c. An expansion for He^ is obtained by 
replacing J 2]c+V on the right-hand side by the Hankel function H 2 i? +P of the second type. As it 
turns out, however, this latter expansion is very poorly convergent for the present problem, 
and it is better to utilize the relation 

He$ (u,v) = Je m (u,v)—iNe m (u,v), 

along with the expansion for the radial Mathieu of the third type, given by (see footnote 6) 

and 

(6) 



Ne 2r+1 (c,u)=(-iy (£f g (-1)* ^g±2 lY k+1 (<x)J k (P)+Y k (<x)J k+1 m, 

where a=ce u /2 and P=ce~ u /2. The functions J k and Y k are Bessel functions of the first and 
second kind, respectively, of order k. The corresponding expansions for the derivatives of the 



6 Tables relating to Mathieu functions (Columbia University Press, New York, N. Y., 1951). 
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Figure 2. 

A , A mplitudc for semicylindrical boss. 

B, Phase for semicylindrical boss. 
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radial function are obtained in a straightforward way by differentiating with respect to u and 
then making use of various recurrence relations for Bessel functions. Therefore, 



pjeW^) -=(f) H (-iy g(-D*Z?^ +y [tanh Ul (2k+ P )J 2k+ ~c sink v t «/« +JF+1 ], 

where the argument of the Bessel functions is c cosh Mi. Correspondingly, 
m meff(fi,u)' 



where 



and 






(7) 

(8) 
(9a) 



r t>Ne 2 T + i(c,u) ~\ 



=(-D r (0" 2 g (-1)* ^ +1 (Y' k+1 (a)J k (p)a+Y' k (a) J k+1 {fi)a 

-Y k+1 (a)J[.i^~Y k (a)J f k + l (m], (9b) 

where the prime over the Bessel functions indicates a derivative with respect to its argument. 
Using the above formulas, the amplitude and phase of the field is computed for several 
cases. The results are normalized so that 

Ae i<p =(H inc '+n aec -)/H , 

where A and <£ are the amplitude and phase of the total field relative to the incident field at 
x=0. 

In figure 2 the functions A and $ are shown plotted along the ground plane in front of 
the ridge at P (i. e., x^>b, y = 0), above the obstacle at Q (i. e,, x=0, y^>a), and to the rear 
of the obstacle at K (i. e., x<^ — b, y=0). In each case, the quantity d indicates the distance 
measured from the elliptical surface to the points P, Q, and R, as illustrated in figure 1. The 
curves in figure 2 refer to the case where pb=fta = 2, so the elliptical boss has degenerated to a 
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Figure 3. 

A, Amplitude of field in front of semielliptical boss. 

B, Phase of field in front of semielliptical boss. 
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Figure 4. 

A, Amplitude of field to rear of semielliptical boss. 

B, Phase of field to rear of semielliptical boss. 
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circular boss whose height is I/t wavelengths. The broken curves in both figures 2, A, and 
2,B, correspond to the trivial situation where the boss is absent. 

It is interesting to note in figure 2 that behind the obstacle, the amplitude of the field is 
reduced in amplitude by about 30 percent and gradually approaches the unobstructed value 
as the observer proceeds away from obstacle. The obstacle imparts an additional phase lag 
of the order of 60° directly behind the obstacle and this value diminishes to about lf)° at two 
wavelengths behind. 

In front of the obstacle, the amplitude and phase of the 4 field oscillates about the unob- 
structed value as the observer proceeds toward the source as a result of the interaction of the 
incoming plane wave with the reflected wave. At a distance greater than a wavelength or 
so, the amplitude varies in the manner of a damped sinusoid with a period of one-half wave- 
length and a magnitude varying inversely as the square root of the distance. In the proximity 
of the obstacle, the structure of the field is quite complex and no simple physical interpretation 
seems possible. It is interesting to note, however, just in front of the obstacle (i. e., x=b, y=0) 
the value of the field is almost twice the value of the incident field. This would be expected 
on the basis of geometrical optics, which predicts a value of exactly two. 

Directly above the obstacle, the field oscillates about the unobstructed value with a 
period approximately equal to a wavelength. The amplitude of the oscillations decreases 
approximately as the inverse square root of the distance. 

To illustrate the effect of changing the cross-seel ion a 1 shape from circular to elliptical, 
values of A and <f> are plotted as of function d both in front (fig. \\) and in the rear (fig. 4) of 
the elliptical boss. Various values of the ellipticity a/6 are shown. Of course, when a/6 
approaches zero the field values approach the unobstructed values. It is interesting to note 
that the qualitative features of the curves are similar, although there is a strong dependence 
of the amplitude and the phase on the ellipticity ratio. In all the calculations, the width of 
the base of the boss (i. e., 2b) has been fixed such that /36 = 2 or 2b is (2/w) of a wavelength. 

The extension of these results to oblique incidence is simple. For example, if the angle 
of the wave front makes an angle 9 with the axis of the elliptical boss, the preceding results 
are applicable if /3 is replaced by sin 6. 

Boulder, Colo., August 23, 1956. 



